Multi-environment testing for G×E interactions and identification of high-yielding, stable, medium-duration pigeonpea genotypes employing AMMI, GGE biplot, and YREM analyses

Pigeonpea [Cajanus cajan (L.) Millspaugh] is a widely grown pulse with high seed protein content that contributes to food and nutritional security in the Indian subcontinent. The majority of pigeonpea varieties cultivated in India are of medium duration (<180 days to maturity), which makes it essential for breeders to focus on the development of stable high-yielding varieties. The diverse agroecological regime in the Indian subcontinent necessitates an efficient multi-environment study by taking into consideration genotype (G) × environment (E) interaction (GEI) that has a significant impact on traits like grain yield (GY) in developing high-yielding and widely adaptable varieties. In the present study, 37 pigeonpea genotypes were evaluated during the 2021 rainy season at ARS Badnapur, ARS Tandur, BAU Ranchi, GKVK Bengaluru, and ICRISAT Patancheru. The GEI was significant on the grain yield (p < 0.01), and hence, genotype + genotype × environment (GGE) and additive main effects and multiplicative interaction (AMMI) biplots along with AMMI stability value (ASV) and yield relative to environmental maximum (YREM) statistics were used to identify stable high-yielding genotypes. The interaction principal component analysis 1 and 2 (IPC1 and IPC2) explained 40.6% and 23.3% variations, respectively. Based on the rankings of genotypes, G37 (ICPL 20205), G35 (ICPL 20203), G8 (ICPL 19404), G17 (ICPL 19415), and G9 (ICPL 19405) were identified as ideal genotypes. Discriminativeness vs. representativeness identified GKVK Bengaluru as an ideal environment for comprehensive evaluation of test genotypes. However, ICPL 19405 was identified as the potentially stable high-yielding genotype for further testing and release across the test environments based on its mean grain yield (1,469.30 kg/ha), least ASV (3.82), and low yield stability index (YSI) of 13.


Introduction
Pigeonpea [Cajanus cajan (L.) Millspaugh] is a grain legume crop of significant economic importance in developing countries in tropical and subtropical regions of the world (Varshney et al., 2012).Globally, pigeonpea occupies an area of 6.03 Mha, producing 5.3 MT with a productivity of 883.4 kg/ha.India contributes up to 65% of the global pigeonpea production with 4.3 MT produced in an area of 5.00 Mha with a productivity of 861.2 kg/ ha (FAOSTAT, 2022).Cooking the dehulled split cotyledons of pigeonpea and serving them as a thick, spicy soup called dal with rice and bread is a common and traditional practice in India.Pigeonpea, containing 21%-25% protein (Saxena et al., 2010), is an ideal complement to cereals for a balanced diet.Beyond its nutritional value, pigeonpea contributes to sustainable agriculture through multiple uses, such as fertilizer (by fixing atmospheric nitrogen and aiding phosphorus release in soil), fuel, fodder, and pharmaceuticals (Mula and Saxena, 2010).Due to rising population demands, there has been an increased need for high pigeonpea seed production, leading to imports valued at nearly 116.57 million USD from African nations (Connect2India, 2021).Currently, the pigeonpea seed production chain constitutes varieties under four major maturity durations, namely, extra-early (90-120 days), early (121-150 days), mid-early (151-165 days), and medium (166-180 days).However, in India, the majority of the pigeonpea cultivated land is dominated by medium-duration maturity varieties.This necessitates the breeder to focus on the development of well-adapted and stable, high-yielding, medium-maturity varieties for the target environments.
The genotype (G) × environment (E) interaction (GEI) plays a significant role in determining the ability of a genotype to thrive in a given environment and ultimately decides its genetic ability to adapt and perform stably.Environment as a whole is a complex of multiple factors that largely include rainfall, temperatures, soil chemistry, soil humidity, disparities in soil type, and biotic stresses like disease and pests, which will together cause GEIs (Oladosu et al., 2016).In this context, extensive evaluation of genotypes across environments is necessary for plant breeding programs that are determined to develop widely adaptable varieties.Multi-environment trials offer the basis for assessing genotypic performance across environments and improving selection accuracy by taking into consideration GEI, which forms the primary factor in this investigation.A multitude of statistical models and tools have been developed to examine the impacts of GEI in mega-environment studies (Eberhart and Russell, 1966).
To facilitate the study of GEI, stability analysis makes use of a variety of statistics.Among them, analysis of variance (ANOVA) dissects the variation into genotype, environment, and genotype × environment effects.Deciphering the stability of the genotypes under study could be conducted by various univariate and multivariate statistical techniques.The best linear unbiased predictors (BLUPs) (Asfaw et al., 2022;Tena and Keneni, 2022;Mwale et al., 2023Mwale et al., , 2023;;Zhang et al., 2023) are estimated for the main and interaction effects of genotype and environment from the combined analysis of variance.It further allows the categorization of the best genotypes according to their environment-specific performance.Many studies have employed combined ANOVA for comparison of data based on multi-environment trials (Lee et al., 2023;Mwale et al., 2023).However, it limits the assessment of genotypic stability, as the model presumes that all genotypic variance and covariance of genotype pairs are the same.Thereby, it forms a restrictive variance-covariance structure of GEI and fails to understand the interaction of each genotype concerning the environments (Hu, 2014).In this direction, multivariate statistics like additive main effects and multiplicative interaction (AMMI) biplots, genotype + genotype × environment (GGE), and yield relative to environmental maximum (YREM) have been developed and employed successfully.
The primary models in GEI analysis are AMMI and GGE biplots (Alizadeh et al., 2017).The AMMI model combines principal component analysis (PCA) and ANOVA into a cohesive approach that may be applied to the analysis of multilocation trials (Zobel et al., 1988;Crossa et al., 1990;Gauch and Zobel, 1996).The GGE biplot uses a two-way table to graphically display a genotype × environment interaction (Yan et al., 2000).This tool is useful for a variety of tasks, including genotype evaluation (the mean performance and stability), test-environmental evaluation, and mega-environment analysis (e.g., "which-won-where" pattern), whereby particular genotypes can be recommended for specific areas (Ding et al., 2007).A unique kind of standardized estimate of a genotype's performance with a negated environment primary effect is called YREM.Additionally, YREM provides an understandable measure of test genotype performance that is independent of genotype attendance (Yan, 1999), which aids in cultivar assessment and serves to identify crossover genotype × environment interaction (Ashwini et al., 2021;Spoorthi et al., 2021).Because each of the aforementioned statistics has merit, they are combined to produce more accurate results.The aim of the present study was therefore to identify the reliable mediumduration pigeonpea genotypes demonstrating stable higher yield using stability models such as AMMI biplot, GGE biplot, BLUP, and YREM.

Genetic material
The experimental material comprised 37 pigeonpea genotypes of medium-maturity duration developed at ICRISAT, Patancheru, and three checks, viz., ICPL 87119, ICPL 8863, and a local check (high-performing variety of that environment).The list of the 37 genotypes used in the study is provided in Supplementary Table 1.

Field trial evaluation
The genotypes along with three checks were evaluated in an alpha-lattice experimental design with three replications during the 2021 rainy across five distinct pigeonpea-growing regions in India, viz., ARS Badnapur (Maharashtra), ARS Tandur (Telangana), BAU Ranchi (Jharkhand), GKVK Bengaluru (Karnataka), and ICRISAT Patancheru (Telangana).The latitude and longitude of the test locations are provided in Supplementary Table 2.Each entry was sown in four rows of 4-m length with a row-to-row spacing of 75 cm and plant-to-plant spacing of 25 cm.The recommended package of practices was followed to raise a healthy crop.Data were recorded on five plants per plot, and observations were taken on five traits, namely, days to 50% flowering (DF) in days, days to maturity (DM) in days, plant height (PH) in cm, grain yield (GY) in kg/ha, and 100 seed weight (HSW) in g.

Statistical analysis
Data on five quantitative traits were subjected to combined ANOVA across five environments to assess the main and interaction effects of genotypes and environments, considering genotypes, environments, replication, and block as random effects.The individual variance of environments was estimated and modeled to error distribution using the residual maximum likelihood (REML) estimator procedure with ASREMLv4.2(Butler et al., 2017).BLUPs were estimated for all main and interaction effects from the combined analysis of variance.
Statistical analysis on AMMI and GGE to create the biplots was carried out in R studio version R 4.1.3(Ahmad et al., 2023).Stability parameters such as AMMI stability value (ASV) and stability index (SI) were estimated based on the AMMI model to assess the relative stability of genotypes using R v4. 1.3 (Chambers, 2008).YREM estimates in the present study were calculated based on BLUP values using Microsoft Excel software.

Combined analysis of variance
Analysis of variance described the fixed and random effects of genotype, environment, and genotype × environment, and the results displayed a significant consequence of environment over the genotypes.Pooled ANOVA revealed that the main and interaction effects of genotype and environment were significant for all the traits in the study including days to 50% flowering (days), days to maturity (days), plant height (cm), 100 seed weight (g), and grain yield (kg/ha) (Table 1).Combined analysis of variance results revealed that the random effects of the environment, genotype, and genotype × environment interaction variance components are statistically significant from zero (p < 0.05).Moving further, the percent contribution from each component toward individual traits dissected through the variance components revealed that the environment contributed the highest variation for traits DF, DM, and PH of approximately 42.01%, 78.34%, and 55.52%, respectively.However, genotype contribution is more for HSW (44.10%), and genotype × environment contributed the highest variation for the GY of approximately 61.42%.With heritability being the key genetic component under selection, among the five traits studied, high heritability was observed for DM (85.91%) followed by DF (84.59%) and HSW (81.02%) (Supplementary Table 3).Correlation between the traits studied showed that DF had a high positive correlation with DM (r = 0.988, p < 0.01) and PH (r = 0.521, p < 0.01).Similarly, DM showed a positive correlation with PH (r = 0.513, p < 0.001) (Supplementary Table 4).

AMMI biplots for assessing stability
Among different AMMI models, widely used AMMI 1 and AMMI 2 were determined in the present study to identify the highly adaptable and stable genotype through multi-environmental trial evaluation considering the genotype main effect and genotype × environment interaction effect.The results of AMMI 1 and AMMI 2 are presented below.

AMMI stability value, stability index, and yield stability index
The genotypes selected based on stability were further analyzed using the AMMI stability value.The ASVs in the present study were estimated using both interactive principal components, IPC1 and IPC2 scores, and it turns out that they significantly contributed to the total genotype × environment variance of grain yield.The ASV estimates are indirectly proportional to the stability of a genotype; i.e., the lower the estimate, the greater the stability.In the present study, the ASV among test genotypes ranged from 1.08 for G11 (ICPL 19407) to 40.72 for G36 (ICPL 20204).Genotypes G9 (3.82), G35 (6.48), G12 (9.26), and G21 (9.27) had the lowest ASVs along with the above trial mean grain yield performance (Table 2).Furthermore, the  Additive main effects and multiplicative interaction (AMMI) biplots of test genotypes for grain yield (kg/ha). SI estimated for test genotypes across environments ranged from very low (13%) for G9 (ICPL 19405) to high (75%) for G34 (ICPL 20202).Among all the best three genotypes with the lowest SI values, G9 (13%) and G35 (24%) had the lowest ASVs as well and performed higher than the trial mean for grain yield.The YSI was determined based on the rank of the mean grain yield of genotypes (RY) and the rank of the AMMI stability value (RASV) for test genotypes, and it was observed that the lowest YSI was observed for G9 (ICPL 19405) of approximately 13 and the highest for G34 (ICPL 20202) of approximately 75.The best genotypes identified based on the mean grain yield, ASVs, and YSI are presented in Table 2.

Genotype × environment mean interaction plot
A combined analysis of variance was used to calculate the genotype × environment means (BLUPs).The mean interaction plot (Figure 2) depicts how the genotype grain yield means were changed in their magnitude across the test environments.Out of 39 genotypes (including checks), 18 were performed above the trial mean wherein G36 (ICPL 20204) performed the highest with 1,441.26kg/ha followed by G14 (ICPL 19411) with 1,419.43 kg/ha and G17 (ICPL 19415) with 1,383.97kg/ha.The 10 best-performing pigeonpea genotypes for mean grain yield are enlisted in Table 3.In contrast, 21 genotypes performed below the trial mean with the lowest of 1,186.47 kg/ha by G18 (ICPL 19416).ICRISAT Patancheru recorded the highest mean value of 1,590.70 kg/ha, and the overall trial mean across the genotypes and environments was 1,318.88kg/ha.

Genotype + genotype × environment biplot
The genotype + genotype × environment interaction biplot for the grain yield trait is represented in Figure 3.Under this, principal component 1 (PC1) contributed 43.13% and principal component 2 (PC2) contributed 20.83% to the total variation.In addition, PC1 and PC2 together contributed approximately 63.96% to the total variation for grain yield.GGE biplot facilitates visual cultivar evaluation and assesses the relative stability of test genotypes.A ranking biplot drawn thus classifies the genotypes based on the environments.Genotypes located near the IPC axis and around the Average environment coordinate (AEC) point are the ideal genotypes.G35 (ICPL 20203), G37 (ICPL 20205), G11 (ICPL 19407), and G9 (ICPL 19405) genotypes had shorter projections from the AEC axis across environments.Meanwhile, G18 (ICPL 19416), G20 (ICPL 19418), and G28 (ICPL 19426) had shorter projections but were placed far away from the AEC arrow (Figure 3A).Environmental vectors were represented on the biplot by perpendicular lines from the origin.The environmental vector of GKVK Bengaluru was the longest as compared to that of other environments.Biplot was divided into mega-environments by lines originating from the biplot origin, based on grain yield data across environments.Test environments were divided into two megaenvironments.BAU Ranchi and ARS Tandur were placed in one mega-environment, while ICRISAT Patancheru, ARS Badnapur, and GKVK Bengaluru were placed in another (within two sectors).Few genotypes occupying the corners of the polygon fall in the sector without any environmental vectors in it, such as ICPL 19416, ICPL 20202, and ICPL 19432.Genotype G36 (ICPL 20204) had a longer projection length from the axis and occupied the vertex of the polygon in the first mega-environment (Figure 3B).In contrast, genotypes G30 (ICPL 19428) and G14 (ICPL 19411) were placed near the environmental vectors of BAU Ranchi and ARS Tandur.The angle between environmental vectors talks about the correlation of results among them.BAU Ranchi and ARS Badnapur vectors formed an obtuse angle (angle > 90°) between them, indicating that the results of the two environmental vectors are negatively correlated.

Estimation of yield relative to environmental maximum model for cultivar assessment
In the present investigation, variation in YREM was evident in each of the environments and was significantly different from others.Variations in YREM and predicted values across environments are depicted in Figure 4.The variation in the estimates was lower in BAU Ranchi, and genotypes showed more variation in ARS Badnapur.The whisker plot of GKVK Bengaluru with lower YREM values occupied a lower part of the graph than other environments.YREM estimates of best pigeonpea genotypes for grain yield across environments are presented in Table 4. Genotypes G14 (ICPL 19411) and G36 (ICPL 20204) both had the highest average YREM value of 0.84 estimated across five environments.Subsequently, genotypes G21 (ICPL 19419), G29 (ICPL 19427), G17 (ICPL 19415), and G12 (ICPL 19408) showed YREM values of 0.77 each, whereas G9 (ICPL 19405) and G30 (ICPL 19428) showed YREM values of 0.75 each.

Discussion
Enhancing the productivity of pigeonpea assumes specific significance in the Indian subcontinent to meet the increasing demand and nutritional security in terms of the daily protein requirement of the predominantly vegetarian population who consume pigeonpea in the form of dry split dal.The genotype × environment interactions, or multi-environment trials, are essential components of crop genetic improvement and breeding programs.

A B FIGURE 4
Box-whisker plots depicting (A) yield relative to environmental maximum (YREM) values of 37 pigeonpea genotypes across test environments and (B) variation in predicted values of genotypes for test environment.

A B FIGURE 3
Genotype + genotype × environment (GGE) biplots for grain yield.According to the combined analysis of variance, the percentage contribution of both environment and genotype × environment interaction effects explained more than 70% of the variation for all traits except HSW, indicating that the individual influence of genotypes did not differ significantly and that there is significant non-crossover interaction that exists (G×E); that is, the genotype mean performance will be affected by the environments.Consequently, the current study uses multi-environmental trials to establish a few stability models to find stable and adaptable genotypes and also to evaluate crossover interactions, such as AMMI, GGE, and YREM.
Additive main effects and multiplicative interaction 1 biplot AMMI is one of the essential models for assessing the impact of genotype × environment interactions on economically important traits like grain yield and its related traits across several environments.Understanding the interplay between genotypes and the relevant environments is made possible by the AMMI model.AMMI 1 is primarily used to discover high potential yield and stability, according to Olivoto et al. (2019).While analyzing AMMI 1, Kılıc ̧(2014) reported that the genotype and environment mean when positioned parallel to the ordinate indicates nearly equal performance.However, a higher yield is displayed by the genotypes positioned on the right side of the biplot's center than by those on the left.These results were found similar to our findings for grain yield.The degree of interaction between genotypes increases, which differs from the origin; the less interactive the genotypes, the closer they are to the origin.In the present study also, genotypes placed closer to the origin were less interactive, suggesting that the genotypes exhibit significant adaptability and possess favorable characteristics for achieving grain yield, whereas genotypes like G36 (ICPL 20204), which was placed farther from the origin, were found to be more interactive, suggesting that these genotypes demonstrated a constrained capacity for adaptation and are better suited for environments characterized by limited conditions.
Within the framework of the AMMI 1 study, it has been noted that genotypes close to the IPC1 axis center region show higher stability with fewer interaction effects.These genotypes exhibit a broad spectrum of adaptation to diverse environmental conditions as a result.In contrast, when the environment and genotype on the IPC axis have congruent polarity, a positive interaction is detected.However, when the environment and genotype have different polarities, an undesirable interaction results.The results obtained from our investigation align with the recorded findings (Kaya et al., 2002;Muniswamy et al., 2021;Reddy et al., 2023;Kona et al., 2024).The results of this study provide empirical support for the implementation of AMMI 2, showing that the AMMI model demonstrates a satisfactory degree of consistency with the gathered data.

Additive main effects and multiplicative interaction 2 biplot
The environmental and genotypic scores of the first two AMMI components are used to create the AMMI 2 biplot.Understanding the role of genotype × environment interaction and the adaptability of the genotypes in the test environments is aided by the IPC1 and IPC2 scores.The first two principal components explained approximately 64.5% of the variation.From the AMMI ANOVA (Supplementary Table 5), the first four IPCs [IPC1 (41.23%),IPC2 (23.20%),IPC3 (20.30%), and IPC4 (15.1%)] explained approximately 99.83% of the variation in genotype × environment TABLE 4 Estimates of yield relative to environmental maximum (YREM) of 10 best pigeonpea genotypes across environments for grain yield.

Genotype
ARS Badnapur ARS Tandur BAU Ranchi GKVK Bengaluru ICRISAT Patancheru Average YREM interactions that were found to be significant.However, several other authors have also reported the predominance of IPC1 and IPC2 interaction components in explaining maximum variation for yield traits (Gauch and Zobel, 1996;Khan et al., 2021;Esan et al., 2023;Kona et al., 2024).
From the AMMI 2 biplot, more general adaptation is expressed by genotypes that are closer to the ordinate axis, and more environment-specific adaptability is expressed by genotypes that are farther from it.AMMI 2 was used in this investigation to identify the specific adaptable genotypes for the respective environments.The results of the current study are consistent with those of Purchase (1997), and also similar findings were reported by Kılıc ̧(2014) and Khan et al. (2021).The environment GKVK Bengaluru exhibits longer vectors, suggesting their greater contribution toward genotype × environment interactions and hence more discriminatory ability of the genotypes evaluated.Genotypes G36 (ICPL 20204) and G37 (ICPL 20205) were found to be more adaptable to the highly discriminating environment of GKVK Bengaluru in the current study.
To quantify and rank genotypes according to their yield stability, the AMMI model does not include a quantitative stability measure.Purchase et al. (2000) developed the ASV measure as a solution to this issue.The genotype stability is indicated by the ASV.Low ASV genotypes are thought to be more stable, whereas high levels indicate less stable genotypes (Hagos and Abay, 2013).Selection is not necessary because a genotype that routinely produces low yields can nonetheless be stable in terms of yield performance (Yan and Tinker, 2006).The most stable genotypes do not always have the best yield performance in certain situations (Oliveira and Godoy, 2006).For this reason, when estimating the YSI, high grain yield is taken into account together with stability.To choose varieties, the YSI combines yield and stability in a variety of settings into a single index.The YSI sums the rank of mean yield across environments with the rank of the ASV of genotypes (Tumuhimbise et al., 2014;Baraki et al., 2014).Genotypes with lower YSI are desirable since they combine high mean yield performance with stability (Tumuhimbise et al., 2014;Baraki et al., 2014;Bose et al., 2014).Based on the YSI, genotypes G9 (ICPL 19405), G12 (ICPL 19408), G14 (ICPL 19411), G21 (ICPL 19419), G17 (ICPL 19415), G29 (ICPL 19427), and G35 (ICPL 20203) were selected as combining high yield performance with stability.

Genotype + genotype × environment biplot
GGE biplot displays genotype + genotype × environment interaction study using site regression analysis, which explains genotype evaluation (stability and adaptability), environment evaluation (representativeness and discriminating power), and mega-environment (which-won-where) evaluation (Yan et al., 2000;Vemula and Parthasarathy, 2023).The first two principal components explained 64% of the total variation.Ideal genotypes are supposedly located near the IPC axis around the AEC point.The relative lengths of projections of the genotypes from AEC are indicative of their relative stability.The shorter the length of the projections of genotypes from AEC, the greater the stability of the genotypes and vice versa (Yan and Kang, 2003).Genotypes G35 (ICPL 20203), G37 (ICPL 20205), G11 (ICPL 19407), and G9 (ICPL 19405) had shorter projections from the AEC axis and were present around the AEC arrow (toward the AEC).These genotypes are widely adaptable, and performance is consistent across the environments.Genotypes G2 (ICPL 19395), G20 (ICPL 19418), and G28 (ICPL 19426) also had shorter projections but were placed in opposite directions of the AEC axis, indicating that they are highly stable and have poor performance of grain yield.Similarly, GGE biplots were also used by researchers to identify stable genotypes in previous studies in pigeonpea (Srivastava et al., 2012;Kumar et al., 2021;Rao et al., 2022).
Each environment is represented as a vector in the biplot by drawing the perpendicular lines from the origin.Longer environmental vectors are indicative of their discriminative ability.GKVK Bengaluru environment showed high discriminating power depicting the variation among the genotypes as high in this environment.Genotypes express themselves, thereby providing a stage to discriminate them based on their performance.Perpendicular lines drawn from the origin divide the biplot into sectors called megaenvironments (Kumar et al., 2021).A mega-environment is made up of a group of environments existing in between the sectors.If any of the genotypes fall in the mega-environment region, then they are considered as specific adaptable (Yan, 1999).BAU Ranchi and ARS Tandur shared similar conditions, while ICRISAT Patancheru, ARS Badnapur, and GKVK Bengaluru environments provided similar growing conditions to the experimental material.The angle formed between the environmental vectors talks about the crossover interaction.When the cosine angle between the environmental vectors is acute (<90°), the results from the two environments are correlated (Yan et al., 2000;Farshadfar et al., 2012), thereby indicating that genotype ranks are not changing across the environments.If the cosine angle between environmental vectors is obtuse (>90°), then genotype ranks are changing across the environments.In the current study, BAU Ranchi and ARS Badnapur environmental vectors formed an obtuse angle; then, the presence of crossover interaction between BAU Ranchi and ARS Badnapur (angle >90°) denoted that the genotype ranks were changing across these two environments.

Yield relative to environmental maximum
To identify crossover genotype × environment interaction and measure the decline in test genotypes' potential grain yield due to crossover genotype × environment interaction, YREM, a straightforward statistic, was employed (Spoorthi et al., 2021).The greater the genotype's YREM value, the smaller the crossover genotype × environment interaction magnitude and the smaller the genotype's potential drop in grain production even in the presence of crossover genotype × environment interaction.The performance of the best genotypes is what it can potentially attain in a particular environment.Therefore, YREM serves as a proxy for the crossover GEI's magnitude.Consequently, the average YREM of a genotype evaluated across environments must equal 1.0 in the absence of crossover genotype × environment interaction.Any deviation of a genotype's YREM from 1.0 is understood as a crossover genotype × environment interaction-related reduction in the genotype's achievable grain yield (Yan, 1999).In the present investigation, YREM values varied across environments, indicating that yield reduction was noticed as the effect of genotype × environment interaction.Among all the test environments, in GKVK Bengaluru, lower YREM values were seen, which is due to the higher yield reduction.A higher YREM value of 0.84 was achieved by G14 (ICPL 19411) and G36 (ICPL 20204) across five environments, indicating that 84% of their grain yield potential can be realized across the test environments, making them relatively stable genotypes.The estimates of YREM were also employed earlier for diagnosing genotype × environment interaction in studies conducted by Yan (2000), Fikere et al. (2008), andAshwini et al. (2021).

Conclusion
Enhancing pigeonpea productivity is crucial for the Indian subcontinent, addressing increasing demand and nutritional needs, particularly for the predominantly vegetarian population.This study utilized multi-environment trials (AMMI, GGE, BLUP, and YREM) to identify medium-maturity genotypes that are adaptable and stable.Combined ANOVA for traits like DF, DM, PH, HSW, and GY indicated significant differences among 37 pigeonpea genotypes, revealing varied performances across environments.Significant genotype × environment interactions highlighted the necessity for in-depth analyses to refine cultivar selection.The environment contributed most to variations in DF, DM, and PH, while genotype × environment had the highest influence on GY.Genotype G36 (ICPL 20204) was identified as widely adaptable with high BLUP values.AMMI and GGE biplots further assessed the impact of genotype × environment interactions.AMMI 1 biplot analysis identified high-yield and stable genotypes, while AMMI 2 provided insights into adaptability using IPC analysis scores.The GGE biplot effectively ranked genotypes and identified ideal test environments.Genotypes G37 (ICPL 20205), G35 (ICPL 20203), G8 (ICPL 19404), G17 (ICPL 19415), and G9 (ICPL 19405) were found to be stable and high-yielding, while environments like GKVK Bengaluru demonstrated high discriminatory power.YREM values were used to quantify crossover genotype × environment interaction effects, with genotypes G14 (ICPL 19411) and G36 (ICPL 20204) showing minimal yield reduction, indicating stability.These comprehensive analyses support the selection of genotypes that combine high-yield performance with stability across diverse environments, aiding in the development of resilient pigeonpea cultivars.

TABLE 1
Pooled analysis of variance of the 37 pigeonpea genotypes across five environments.
* and ** are significant at the probability of 0.05 and 0.01, respectively.The percentage contribution of main and interaction effects of environment and genotype is in parentheses.Numbers highlighted in bold are the highest contributions.DF, days to 50% flowering; DM, days to maturity; PH, plant height; HSW, 100 seed weight; GY, grain yield.

TABLE 2
Estimates of mean grain yield and additive main effects and multiplicative interaction (AMMI) model-based parameters to assess the stability of the pigeonpea genotypes for grain yield.

TABLE 2 Continued
Genotypes with bold letters indicate the best genotypes with the above trial mean grain yield having ASV less than 10 and lowest YSI.RY, rank of the test genotype based on mean grain yield; ASV, AMMI stability value; RASV, rank of the test genotype based on ASV; SI, stability index; YSI, yield stability index.

TABLE 3
Estimates of genotype, environment, and genotype × environment BLUP values of the 10 best genotypes for grain yield in pigeonpea.